

EFFECT OF LIMB DARKENING ON EARTH RADIATION 
INCIDENT ON A SPHERICAL SATELLITE 

by S. Katzoff and G, Louis Smith 

Langley Research Center 
Hampton f Va. 2S665 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • DECEMBER 1974 






1. Report No. 2. Government Accession No. 

NASA TN D-7807 


4. Title and Subtitle 

EFFECT OF LIMB DARKENING ON EARTH RADIATION 
INCIDENT ON A SPHERICAL SATELLITE 


7. Autbor(s) 

Se Katzoff and G, Louis Smith 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


15, Supplementary Notes 


3. Recipient's Catalog No. 


5. Report Date 

December 1974 


6* Performing Organization Code 


8. Performing Organization Report No. 

L-9698 



11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Note 


14. Sponsoring Agency Code 


16. Abstract 


The thermal radiation from the Earth incident on a spherical satellite depends on 
the angular distribution of Earth -emitted radiation. (The deviation from a Lambertian 
distribution is commonly referred to as limb darkening.) An analysis is presented of 
this dependency, and calculated results are given, based on a published limb -darkening 
curve for the Earth. The curve was determined from Tiros data, and is a statistical 
average over the entire globe between ±75*^ latitude. The computed effect of limb dark- 
ening was 1.8 percent at 900 km altitude, 2.5 percent at 500 km altitude, and 3.0 percent 
at 300 km altitude. Below 300 km, it increased rapidly with decreasing altitude. 

Discussion is included of various other problems inherent in the use of orbiting 
spheres and stabilized flat plates to measure the heat radiated from the Earth. 


17. Key Words ^Suggested by Author Is}) 


1 18. Distribution Statement 


Energy budget of Earth 
Radiation 
Satellite heating 
Spherical satellites 
Earth limb darkening 
Earth thermal radiation 


Unclassified - Unlimited 

STAR Category 31 

19. Security Dassif. (of this report) 

20, Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 


23 

$3.25 


For sale by the National Technical Information Service, Springfield. Virginia 22151 





















EFFECT OF LIMB DARKENING ON EARTH RADIATION 
INCIDENT ON A SPHERICAL SATELLITE 

By S, Katzoff and G. Louis Smith 
Langley Research Center 

SUMMARY 

The thermal radiation from the Earth incident on a spherical satellite depends on 
the angular distribution of Earth-emitted radiation. (The deviation from a Lambertian 
distribution is commonly referred to as limb darkening.) An analysis is presented of 
this dependency, and calculated results are given, based on a published limb- darkening 
curve for the Earth, The curve was determined from Tiros data, and is a statistical 
average over the entire globe between ±75^ latitude. The computed effect of limb dark- 
ening was 1.8 percent at 900 km altitude, 2.5 percent at 500 km altitude, and 3.0 percent 
at 300 km altitude. Below 300 km, it increased rapidly with decreasing altitude. 

Discussion is included of various other problems inherent in the use of orbiting 
spheres and stabilized flat plates to measure the heat radiated from the Earth, 

INTRODUCTION 

This paper discusses the following problem: Given the intensity (power per unit 
area) of the thermal radiation emitted from the Earth, to what extent is the corre- 
sponding irradiation of a spherical satellite a function of the directional distribution of 
this emitted radiation? The problem arises because it has been proposed to study the 
energy budget of the Earth through measurements of the temperatures within a fleet of 
such orbiting spheres; for this purpose, it is obviously necessary that the temperatures 
be essentially independent of the directional distribution of the Earth-emitted radiation, 
or else it is necessary to account for such a dependency. 

The analysis is made for the cases in which the directional distributions are given 
as (1) a power series in the zenith angle, and (2) a power series in the cosine of the 
zenith angle. In order to show the order of magnitude of the effect for a practical case, 
some calculated results are given in which experimental satellite data for the directional 
distribution were used. 



The following section provides some further background concerning the special 
purpose of the thermal balance measurements, the desired accuracy, and some possible 
approaches, with their inherent difficulties. It should also help to put the subject of the 
present analysis in perspective as one among several troublesome basic questions 
regarding the use of the sphere as an integrating radiation sensor in attempting to meas- 
ure the Earth’s thermal balance to high accuracy. 

PERSPECTIVE OF THE PROBLEM 

Measuring and understanding the Earth’s energy budget have been among the basic 
problems in meteorology. The recent development of popular concern with the environ- 
ment, however, has made more general the recognition that gradual modification of the 
atmosphere and of land and water areas, resulting both from natural causes and, espe- 
cially, from man's activities, may eventually produce significant, or even profound, 
changes in the Earth’s energy budget and, hence, in the Earth's climate. With the devel- 
opment of space technology, meteorologists recognized the feasibility of using satellites 
for measuring and even continuously monitoring the Earth's energy budget and also 
regional energy budgets; and they have accordingly pushed the development of the neces- 
sary technology. 

High accuracy is necessary if slow modifications and long-term trends are to be 
recognized - for example, for albedo measurements, the accuracy requirements as 
estimated by various groups range from 0.1 percent to 1 percent. A fleet of satellites 
capable of jointly observing all areas of the Earth nearly simultaneously would obviously 
be needed, but the nature of their measurements remains under discussion, Horizon-to- 
horizon scanning of the Earth, with appropriate spectral resolution, spatial resolution, 
and radiometric accuracies, would seem to provide the most thorough and unequivocally 
interpretable information. Indeed, our present estimates of the Earth's albedo are sup- 
ported by just such data. However, the dubiousness of achieving long-term reliability 
of the equipment and instruments, and the practical difficulties in working up huge quan- 
tities of data militate against such an approach for a continuous long-term project. 

Integrating Radiometer -Type Sensors 

Other approaches, presumably simpler and more practical, involve integrating 
radiometer -type sensors that respond to the total radiation received from the spherical 
segment of the Earth below them (horizon-to -horizon in all directions). A group of 
two or three such sensors on each satellite, differently coated so that they have 
different absorptances for sunlight (direct or Earth-reflected) and also have different 


2 



absorptances for Earth -emitted thermal radiation, should provide a set of simultaneous 
equations from which the separate intensities of these two wavelength bands can be 
extracted. Two such sensor concepts will be briefly discussed in the following 
paragraphs. 

The flat-plate sensor .- Consider an imaginary spherical surface concentric with 
and surrounding the Earth, as in sketches (a) and (b). All the radiation leaving the 
Earth, either as reflected sunlight or as Earth emission, must pass outward through 
this spherical surface; this must be true regardless of the directional distribution of the 
emitted radiation - for example, whether the radiation is purely radial, as in sketch (a), 
or is diffuse, as in sketch (b). A small flat-plate sensor p in a circular orbit, sta- 
bilized such that it is always normal to the radius, may be considered as an element of 



Sketch (a) - Radial emission Sketch (b) - Diffuse emission 

this imaginary surface, intercepting the radiation passing outward at that instant at that 
location. An adequate number of such sensors, all in circular orbits of the same radius, 
could thus represent the surrounding imaginary sphere and could, accordingly, adequately 
sample the outgoing radiation from the Earth, regardless of its directionality. The flat- 
plate concept thus seems to be the most obvious approach to the integrating radiometer - 
type sensor mentioned in the preceding paragraph. 

Achieving long-term reliability for the stabilization mechanisms of a fleet of such 
sensors may seem rather formidable, but it is not prohibitive. A rather more basic 
difficulty is caused by the fact that the Earth-reflected solar radiation has a somewhat 
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variable spectral content that, in any case, differs from that of the direct solar radia- 
tion; so that unless the spectral absorptance of the coating is uniform over the entire 
wavelength range of the solar spectrum {which is unlikely), a certain degree of inaccu- 
racy is introduced by assuming that the coating has a unique solar absorptance, appli- 
cable to both direct and reflected sunlight. A similar problem, although probably less 
serious, exists for the infrared, since the Earth-emitted thermal radiation also has a 
variable spectral content. 

A further basic difficulty would result from nonuniform directional absorptance of 
the sensor coating; While it is true that the total radiation striking the sensor is the 
same in sketches (a) and (b), if absorptance decreases with increasing incidence angle, 
less of the radiation will be absorbed in sketch (b) than in sketch (a), so the readings 
will be different. Absorptance does, in fact, generally decrease with increasing inci- 
dence angle, especially for the longer wavelengths. Hence, accurate interpretation of 
the readings would require not only that the directional absorptances of the sensor 
coatings be known over the entire spectral range of interest, but also that the directional- 
spectral distribution of the radiation striking the sensors also be approximately known 
or estimable. This latter requirement is not the same as requiring that the results be 
known beforehand; but it does emphasize the need for special analytical studies to deter- 
mine the best means of taking this matter into account in interpreting the data, and to 
estimate the corresponding range of uncertainty in the answer. 

The spherical sensor .- Instead of a group of flat plates, one might use a group of 
thin-walled hollow spheres, well separated from each other. The basic equations for 
reducing the data from groups of such spheres are given in reference 1. The approach 
is attractive because a uniformly coated sphere does not have to be stabilized, since it 
always presents the same aspect to the Earth regardless of its orientation. Further- 
more, the directional-absorptance question that was raised for the flat plate seems less 
important for the sphere, since the sphere presents the same aspect to radiation received 
from the nadir as to radiation received from the horizon. However, since the angles of 
incidence of any radiation striking the sphere range from 0*^ to 90*^, one must again have 
at least an approximate a priori knowledge of the spectral content of the radiation 
received from the different parts of the Earth's segment beneath it. (Availability of 
experimental data on the directional-spectral absorptance characteristics of the different 
coatings used on the spheres is assumed.) In any case, the problems introduced by the 
nonuniform spectral content of the Earth-emitted radiation and of the Earth-reflected 
radiation, previously discussed for the flat plate, also exist for the sphere. 

Still another complication is introduced if the coating is degraded by solar ultra- 
violet and particle radiation, with corresponding changes in absorptance characteristics. 
Reference 1 indicates that absorptance values may be extracted from the flight data, so 
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the degradation might seem to be only an acceptable annoyance; however, unless the sur- 
face of the sphere is uniformly exposed to the degrading radiations, the charges in absorp- 
tance characteristics will not be uniform over the surface. Uniform exposure might be 
provided by using magnetic coils in the satellite to interact with the Earth's magnetic 
field and suitably rotate the satellite, but the basic simplicity of the satellite would 
thereby be compromised. 

A basic source of uncertainty in interpretating the data is that the irradiation by 
the Earth depends on the directional distribution of the radiation leaving the Earth. It 
will be remembered that for the flat-plate sensor the amount of radiation received at 
the flat plate from the Earth is independent of this directional distribution - the flat 
plate is merely a sample area of the imaginary surrounding sphere (sketches (a) and 
(b)) and intercepts all the radiation that would be passing outward through that area, 
regardless of the direction from which the radiation comes. For the spherical sensor, 
a similar statement caivnot be made. If the flat-plate sensor is replaced with, a spher- 
ical sensor in sketches (a) and (b), it would no longer be true that the radiation reaching 
the sensor is the same in both sketches. In fact, in a theoretical limiting case, when the 
sphere is just above the Earth and the directional distribution of radiation in sketch (b) 
is Lambertian, the sphere would receive twice as much radiation as it would for the 
purely radial radiation of sketch (a). 

Fortunately, this 2:1 uncertainty is not realistic, since (1) the sphere would not be 
very close to the Earth, because the requirements for a long lifetime and a broad over- 
view require that the satellite be at a fairly high altitude (of the order of, say, 1000 km), 
and (2) the purely radial distribution in sketch (a) never occurs, and we are, in fact, 
concerned only with moderate deviations from the Lambertian distribution. Both Earth- 
emitted thermal radiation and, especially, Earth- reflected solar radiation show devia- 
tions from Lambert' s law. 

The present study concerns specifically the irradiation of a sphere by the Earth- 
emitted radiation, for which the deviation is often referred to phenomenologically as 
limb darkening. As in other typical cases of limb darkening, it is characterized not 
only by reduced infrared brightness of the horizon (compared with the nadir) but also 
by relatively more of the longer wavelengths in the radiation from the horizon; and it 
exists because the radiation from the horizon traverses a long path through the cold 
upper troposphere and lower stratosphere, whereas the radiation from the nadir tra- 
verses a shorter path through these cold air layers and, accordingly, remains more 
nearly like that originally emitted from the warm Earth’s surface and lower air layers. 

Reference 2 called special attention to the effect of limb darkening on the input to 
a spherical sensor; and, since then, the Tiros series of satellites has provided data 
from which the limb darkening of the Earth- atmosphere system has been deduced 
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(refs. 3 and 4). The data from references 3 and 4 have been used in this study to calcu- 
late the limb-darkening effect on the irradiation of a spherical satellite. In addition, 
as mentioned previously, curves and formulas are given herein that may be applied for 
any specified deviation from the Lambertian radiation law. In the application of these 
results, it may be desirable to allow for the fact that the wavelengths of the radiation 
from the limb are larger than those from the nadir; because if the directional absorp- 
tance characteristics of the sphere coatings are not constant throughout the entire 
infrared range of interest, some further correction for limb darkening would be needed. 


SYMBOLS 

A cross-sectional area of spherical satellite, 

, , 0" sin » cos » ,, „„„ 

sin^^ 

a coefficient in power series for F{(f)) (see eq. (8)) 


D 


factor measuring effect of limb darkening on Earth’s thermal radiation 


incident on spherical satellite, 


I(k) 

An(k) 


F(0) 


limb -darkening function 


"n 


coefficient in cosine power series for F{4>) (see eq. (A2)) 


h satellite altitude, km 

I(k) = F( 0) stn 0 cos 0 

do . f ' ^ O 


k = sin 9 




sin^ <p 


intensity of radiation (per steradian) from a surface point (per unit 
projected area), W/m^-sr (eq. (1)) 


Q 


irradiation of spherical satellite by Earth, W (eq. (3)) 


R 


radius of Earth, km 
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s 


satellite 


T temperature of Earth, K 

a nadir angle of point on Earth’s surface from satellite, radians (see fig. 1) 

e total hemispherical emissivity of Earth (considered to be uniform over area 

seen from satellite) 

9 angle between nadir and tangent to Earth’s surface from satellite (nadir 

angle of horizon, see fig. 1), radians 

a Stefan -Boltzmann constant, 5.6696 X 10"^ w/m^-K"^ 

(p angle between zenith and line from Earth’s surface to satellite (zenith 

angle, see fig. 1), radians or deg 

ANALYSIS AND DISCUSSION 

Basic Equations 

The Earth-atmosphere system will be approximated by a sphere of radius R. The 
radiation leaving any point on the sphere will be assumed to be a function of the zenith 
angle, and independent of azimuth. The intensity of radiation (per steradian) from a 
surface point (per unit projected area) may be written as (see ref. 5) 

1 = - eT^ F(d)) (1) 

7 [ 

where F{4>) is the limb-darkening function and indicates the angular distribution of the 
observed surface brightness. It is a constant, F((^>) = 1, for a Lambertian radiator; and 
in any case, it satisfies the equation 

2 \ F(0) sin (p cos 1 (2) 

•-'0 

Limb- darkening curves for the Earth are given in references 3, 4, and 6. These curves 
are proportional to F(i^), but are normalized such that they are unity at 0=0. 
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Now consider a spherical satellite S at height h above the Earth, as shown in 
figure 1. The radiation from the Earth incident on the satellite is 

Q = A r idf2 = A r i(27f sin a da) 

- Q 

= 2AeeT^ \ F(<^) sin a da (3) 

where dfi is the solid angle subtended at the satellite by an element of the Earth's 
surface and A is the cross-sectional area of the satellite. It is assumed that the 
Earth’s emittance e and temperature T are constant over the view field. The law 
of sines applied to triangle OSP gives 

sin g _ sin 4> 

R R + h 


Using this relation to eliminate a from the integrand in equation (3) gives 

^ sin^ 0 

where 


k = ~~ = sin 0 
R+h 


Denote the integral in equation (4) by I{k) 
nTi/2 (p cos 4> 


I(k) = 


0 \/i^k^ sin^ (p 


d(p 


(4) 


(5) 


( 6 ) 


It is seen that I(k) is a function of k and a function of the limb-darkening curve 
F((p). The influence of limb darkening on Earth irradiation of a spherical satellite is 
contained in I(k). Comparison of equations (2) and (6) shows that 

1 ( 0 ) 4 

Thus, at large distances from Earth, the incident radiation is independent of the limb- 
darkening curve, as would be expected, since at large distances the Earth would approxi- 
mate a point source. 
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Polynomial Limb -Darkening Function 
In reference 3, the limb-darkening curve F’(<]i)) is given as 

F'(<^>) = 1 - 1,116 X 10-4 (f) - 5.873 X 10-5 ^2 

+ 1.387 X 10-6 ^3 _ 1.523 x lO"® 

where (p is expressed in degrees and where the prime indicates that F has been 
normalized such that it is unity at p= 0. 

Because of the second term in this expansion, a plot of F’(4*) against p> has a 
finite slope at 0=0. The existence of such a finite slope seems rather unlikely; for 
example, if F' is plotted against p as the line SP (fig. 1) swings in a plane through 
the zenith, there would be a finite discontinuity in the rate of change of F’ as the line 
passes through the zenith (0 = 0). However, the given equation represents a best fit to 
the data, and altering it in order to eliminate the linear term would hardly be justified. 

In any case, the finite slope at 0 = 0 is quite small. 

The expression for F’(0), when normalized according to equation (2), becomes 

F(0) = 1,0949 - 0.007000- 0.211302 + 0.285603 - 0.179704 (7) 


where p is expressed in radians. 

Following the form of equation (7), denote 


F(0) =^an0*" 
n 

Then equation (6) reduces to 

I(k) = YanAn(k) 
n 


where 



( 8 ) 


(9) 


(10) 
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For the cases k = 0 and k = 1, recursion relations may be readily developed for 
evaluating this integral. Thus, for k = 0, two successive integrations by parts give the 
relation 

An(0) = i(f [ - An-2(0) (11) 

which, together with the easily derived 
Ao(0) = 1 1 
A,(0) = I ^ 

provide an easy recursive determination of all An(0). 

Similarly, for k = 1, two successive integrations by parts give the relation 

An(l) = n(^|) - n(n - l)An_ 2 (l) 

which, together with the easily derived 


( 12 ) 


(13) 


Ao(l) = 1 
Ai(l) = 1 


(14) 


provide an easy recursive determination of all Ai^(l). Table I lists An(0) and An(l) 
for n = 0 to 9. 

The case k = 1 is a limiting case for zero altitude. In reality, radiation from the 
atmosphere would render this analysis inaccurate below an altitude of 150 km. 

For 0 < k < 1, the values of An(k) have been evaluated numerically and are 
shown in figure 2 for values of n from 0 to 9. By use of figure 2, I(k) may be evalu- 
ated for any limb -darkening expression of the form of equation (8). It is seen that the 
values of An(k) vary weakly with k for small k, that is, for large altitudes, and 
vary rapidly for k near 1, that is, for low altitudes. 

An alternative method of evaluating l(k) by expanding F{<p) as a cosine power 
series is presented in the appendix. 
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Results Based on Tiros Data 

The limb -darkening influence bn I(k) corresponding to equation (7) has been 
computed as a function of k and is shown as the solid curve in figure 3, Comparison 
with the upper curve of figure 2(a) shows that I(k) and AQ(k) are nearly the same for 
small k (large satellite altitude) but that the difference becomes appreciable at large 
k (small satellite altitudes). In the absence of limb darkening, I(k) would be identi- 
cally AQ(k). A suitable measure of the effect of limb darkening on heat incident on a 
spherical satellite is the factor 


D = 


I(k) 

A^(k) 


This factor D is plotted as a function of k in figure 4, and as a function of altitude in 
figure 5 (solid curve in both figures). It is seen that D varies rapidly with altitude 
below 300 km although, of course, the lower altitudes would not be practical because of 
air drag. Between 500 and 900 km, the variation is only from 0.976 to 0.984, less than 
1 percent. In order to reduce the effect of limb darkening to less than 1 percent, that 
is, to get a value of D above 0.99, it is necessary to go to an altitude above 2000 km. 

The foregoing analysis has shown the main effect of limb darkening on incident 
radiation from the Earth to a spherical satellite. Some remarks and additional calcula- 
tions concerning the uncertainty of the results are given in the following section. 


Uncertainty of the Limb-Darkening Function 

Results which were presented in figures 3 to 5 were computed by using the limb- 
darkening curve of reference 3, which is a globally averaged curve. Reference 4 pre- 
sents limb -darkening curves for various regions of the Earth: arctic, antarctic, desert, 
tropics, and the southern hemisphere. These curves were normalized in accordance 
with equation (2) and the corresponding incident -heat factors I(k) were computed. 
Although these curves show appreciable differences, the corresponding incident -heat 
factors I(k) for the most part show only very slight differences, and accordingly only 
one I(k) curve is shown in figure 3 for reference 4. For k =0.5, the values of I(k) 
computed for the various curves are equal, to three decimal places. For 0.5 < k < 0,95 
(300 to 6000 km), the values of l{k) differ by no more than 0.007, but for k = 1, the 
values of I(k) differ by as much as 0.03. 

In reference 3, the limb -darkening curves were terminated st ({) = 78.75°. For 
the present analysis, the curve fit (eq. (7)) was needed and was used out to 4> = 90°. 
Also, some of the curves from reference 4 were extrapolated to ^ = 90®. The portion 
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of the curve from 4> = 75° to 90° is relatively more important at low satellite altitudes 
than at high altitudes, because the relative contribution from the areas near the horizon 
is greater. Considerations of the previous paragraph indicate that a moderate inaccu- 
racy in extrapolation is not critical for k < 0.95. 

In figures 3 to 5, it is seen that the incident heat as computed with the extrapolated 
curves of reference 3, which is a global average, is less than the incident heat as com- 
puted with the limb -darkening curves of reference 4. The limb -darkening curve of ref- 
erence 3 lies in the midst of the curves of reference 4, for 0 < 70*^, but drops below 
them for <j> < 80°. This points out the need for measurements of the limb darkening in 
the range 70° < 0 < 90°. 

It was assumed in this study that upwelling radiation is axisymmetric about the 
vertical. When the satellite is at a fairly high altitude, however, the radiation tempera- 
ture of the viewed area may be very nonuntform. For example, if its altitude exceeds 
500 km, the satellite can receive radiation from both the Torrid Zone and a Frigid Zone 
simultaneously. This radial and azimuthal nonuniformity would complicate the interpre- 
tation of data even if the Earth were a Lambertian radiator; but limb darkening adds a 
further complication. Finally, the analysis assumed that the satellite was at a clearly 
definable altitude h above a radiating spherical surface. In reality, radiation is also 
received from clouds and the atmosphere itself, so that the location of the effective 
radiating surface is not clearly definable, especially at the horizons. 

Thus, the present study involves a number of unanalyzed simplifications, investi- 
gation of which could entail considerable effort. Nevertheless, the present analysis 
shows the magnitude and the main characteristics of the limb-darkening effect on Earth 
irradiation of a satellite. 


CONCLUSIONS 

A simple analysis has been made of the effect of limb darkening on Earth irradia- 
tion of a spherical satellite. For this study, an estimated average limb darkening based 
on Tiros data was used. It was concluded that; 

1. Limb darkening must be considered for accurate measurements of Earth - 
emitted heat incident on a spherical satellite, 

2. The computed effect of limb darkening was 1.8 percent at 900 km altitude, 

2.5 percent at 500 km altitude, and 3.0 percent at 300 km altitude. Below 300 km, it 
increased rapidly with decreasing altitude. 

3. In order to reduce the effect of limb darkening to less than 1 percent, it is 
necessary to go above 2000 km altitude. 

4. Measurements of limb darkening in the range 70° < (f) < 90° are needed. 
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The basic assumption that the satellite views a well-defined spherical cap, radi- 
ating uniformly according to a known limb-darkening curve, is recognized as somewhat 
unrealistic. Accordingly, although the present study is adequate for estimating the 
importance of limb darkening, more work is clearly needed to establish a basis for 
relating the Earth thermal radiation to what is received and absorbed by an orbiting 
sphere. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va,, November 8, 1974. 



APPENDIX 


ALTERNATIVE EVALUATION OF I(k) 


In the evaluation of the effect of limb darkening on the Earth's heat incident on a 
spherical satellite, the expression 


Kk) = F(0) sin <t> cos 0 

\ji - fe2 sin2 0 


(Al) 


appears. In the main text, T(<p) was expressed as a power series in consistent 
with Lienesch and Wark (ref. 3). An alternative procedure is to express F(<^) as a 
cosine power series, in which case the integrals required to express the heat incident on 
a satellite can be evaluated in closed form. Thus, let 


F(</>) = ^ Gn cos*i <p 


(A2) 


n=0 


whence 


I(k) = ^ Gn Ln(k) 


(A3) 


where 


n :=0 


^ ,, , sin d> cos'^'*’^ d> , 


^ \/l - k^ sin2 0 

The law of sines (see p. 8) is used to rewrite the integral in terms of a: 


(A4) 






(k^ - sin^ sin o; dof 


A change of variables is made, 


X = cos a 


so that 


" j^n+2 Jb 


n/2 


dx 


(A5) 
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APPENDIX - Concluded 


where 

b = cos 0 = \jl - 

Integration of equation (A5) by parts gives 


(n + 1)L„ = X . . b2) ^ dx 


(A6) 


(A7) ' 


Equation (A5) may be integrated for cases n = 0 and n = 1 to give 


Lo(k) 

LjL(k) 


1 

1 + b 



k + b^ 





(A8) 


The may be computed recursively by using equation (A7), with equations (A8) pro- 
viding the starting values. 


Because F(0) is defined over a hemisphere, another potentially useful form for 
F{(p) is an expansion in terms of Legendre polynomials in cos <p. The Legendre poly- 
nomial may be interchanged with the power series, for example, through use of a table 
as in reference 7, permitting the results of this appendix to be used. 
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TABLE I. VALUES OF A^(0) AND A^(l) FOR n = 0 to 9 


n 

A„(0) 

And) 

0 

0.5000 

1.0000 

1 

,3927 

1.0000 

2 

,3668 

1.1416 

3 

.3799 

1.4022 

4 

.4214 

1.8040 

5 

.4913 

2.3964 

6 

.5947 

3.2576 

7 

.7404 

4.5020 

8 

.9410 

6.3461 

9 

1.2277 

9.4382 
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(b) n = 3, 4, 5, and 6. 
Figure 2.- Continued. 
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(c) n = 7, 8, and 9. 
Figure 2 .- Concluded. 
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Figure 3.- Incident-heat factor I for spherical 
satellite as a function of k. 



Figure 4.- Factor D measuring effect of limb 
darkening on Earth radiation incident on 
spherical satellite as a function of k. 




I 


Limb darkening from ref. 3 (Global average) 


Limb darkening from ref. 4 


Arctic 


Satellite altitude, km 

Figure 5.- Factor D, D = i the effect of limb darkening on Earth radiation incident 

Aq 

on a spherical satellite as a function of altitude. 




